function histmult(x1,x2,binedge,tit); 
% Function plots the histogram of x1 and x2 when not 
% of the same dimension 
% by counting the elements of x contained in each bin 
% whose upper edge is determined by binegde 
% into.title allows to define a title string 
% into.scale =1 divides by the number of rows in x for plotting
[N1,BIN1]=histc(x1,binedge);
figure;
subplot(2,1,1); 
bar(binedge,N1,'histc','r');
subplot(2,1,2); 
[N2,BIN2]=histc(x2,binedge); 
bar(binedge,N2,'histc','b');
if nargin > 3; 
    title(tit); 
end; 
